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Today ion traps are among the most promising physical systems for constructing a quantum 
device harnessing the computing power inherent in the laws of quantum physics 1 ' 2 . The stan- 
dard circuit model of quantum computing requires a universal set of quantum logic gates for 
the implementation of arbitrary quantum operations. As in classical models of computation, 

On 

quantum error correction techniques 3 ' 4 enable rectification of small imperfections in gate 
operations, thus allowing for perfect computation in the presence of noise. For fault-tolerant 
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computation 5 , it is commonly believed that error thresholds ranging between 10~ 4 and 10 2 
^ ■ 6-8 will be required depending on the noise model and the computational overhead for re- 
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alizing the quantum gates. Up to now, all experimental implementations have fallen short 
of these requirements. Here, we report on a M0lmer-S0rensen 910 type gate operation en- 
tangling ions with a fidelity of 99.3(1)% which together with single-qubit operations forms a 
universal set of quantum gates. The gate operation is performed on a pair of qubits encoded 
in two trapped calcium ions using a single amplitude-modulated laser beam interacting with 
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both ions at the same time. A robust gate operation, mapping separable states onto maxi- 
mally entangled states is achieved by adiabatically switching the laser-ion coupling on and 
off. We analyse the performance of a single gate and concatenations of up to 21 gate opera- 
tions. The gate mechanism holds great promise not only for two-qubit but also for multi-qubit 
operations. 

For ion traps, all building blocks necessary for the construction of a universal quantum 
computer 1 have been demonstrated over the last decade. Currently, the most important challenges 
consist in scaling up the present systems to higher number of qubits and in raising the fidelity 
of gate operations up to the point where quantum error correction techniques can be success- 
fully applied. While single-qubit gates are easily performed with high quality, the realisation of 
high-fidelity entangling two-qubit gates 1116 is much more demanding since the inter-ion distance 
is orders of magnitude bigger than the characteristic length scale of any state-dependent ion-ion 
interaction. Apart from quantum gates of the Cirac-Zoller type 2 ' 12 where a laser couples a single 
qubit with a vibrational mode of the ion string at a time, most other gate realizations entangling ions 
have relied on collective interactions of the qubits with the laser control fields 11 ' 13 ~ 15 . These gate 
operations entangle transiently the collective pseudo-spin of the qubits with the vibrational mode 
and produce either a conditional phase shift 17 or a collective spin flip 9 ' 10 - 18 of the qubits. While the 
highest fidelity F = 97% reported to date 13 has been achieved with a conditional phase gate acting 
on a pair of hyperfine qubits in 9 Be + , spin flip gates have been limited so far to F ps 85% 1114 . 
All of these experiments have used qubits encoded in hyperfine or Zeeman ground states and a Ra- 
man transition mediated by an electric-dipole transition for coupling the qubits. While spontaneous 
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scattering from the mediating short-lived levels degrades the gate fidelity due to the limited amount 
of laser power available in current experiments, this source of decoherence does not occur for op- 
tical qubits, i.e. qubits encoded in a ground state and a metastable electronic state of an ion. In the 
experiment presented in this paper where the qubit is comprised of the states \S) = Si/ 2 (m = 1/2) 
and \D) = D 5 / 2 ( m = 3/2) of the isotope 40 Ca + , spontaneous decay of the metastable state reduces 
the gate fidelity by less than 5 • 1CT 5 . 

A M0lmer-S0rensen gate inducing collective spin flips is achieved with a bichromatic laser 
field with frequencies u± = u ± 6, with u the qubit transition frequency and 5 close to the 
vibrational mode frequency v (see Fig. [T]). For optical qubits, the bichromatic field can be a 
pair of co-propagating lasers which is equivalent to a single laser beam resonant with the qubit 
transition and amplitude-modulated with frequency 5. For a gate mediated by the axial centre- 
of-mass (COM) mode, the Hamiltonian describing the laser-qubit interaction is given by H = 
me-^S+ie-W+V + e«(^0) e ^(«-»'+«te»') +k ^ Rere g. = + af^ , j e {+,-,x,y,z}, 

denotes a collective atomic operator, and a+\S)i = \D)i. The operators a, a* annihilate and cre- 
ate phonons of the COM mode with Lamb-Dicke factor 77. The optical phase of the laser field is 
labeled (ft, and the phase £ accounts for a time difference between the start of the gate operation 
and the maximum of the amplitude modulation of the laser beam. In the Lamb-Dicke regime, and 
for = 0, the gate operation is very well described by the propagator 19 

U(t) = e- jF W & J D(a(t)^)exp(-z(At + X sm(z/-5)t)^). (1) 

Here, the operator to the right describes collective spin flips induced by the operator S y ^ = 

S y cosip + 5 z sin^, V = x cos C> and A w tf9?j{y - 5), X « tf^jiy - 5) 2 . With a(t) = 
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a (e l ^~ s ^ — 1), the displacement operator D{f3) = e /3at " /3 * a accounting for the transient entangle- 
ment between the qubits and the harmonic oscillator becomes equal to the identity after the gate 
time T ga te = • The operator e - iF ^ Sx with F(t) = (2Q/5)(sm(5t + 0) - sin0) describes 
fast non-resonant excitations of the carrier transition that occur in the limit of short gates when 
f2 <C 5 no longer strictly holds. Non-resonant excitations are suppressed by intensity- shaping the 
laser pulse so that the Rabi frequency Q(t) is switched on and off smoothly. Moreover, adiabatic 
switching makes the collective spin flip operator independent of ( as S y ^ — > S y for VL — > 0. 
To achieve adiabatic following, it turns out to be sufficient to switch on the laser within 2.5 trap 
cycles. In order to realize an entangling gate of duration r gate described by the unitary operator 
Ugate = exp(— i^Sy), the laser intensity needs to be set such that r]Q, \8 — H/4. 

Two 40 Ca + ions are confined in a linear trap 20 with axial and radial frequencies of Vaxiai/^-n = 
1.23 MHz and Vradiai/^K = 4 MHz, respectively. After Doppler cooling and frequency-resolved 
optical pumping 21 , the two axial modes are cooled close to the motional ground state (h CO m, 
hstretch < 0.05(5)). Both ions are now initialized to \SS) with a probability of more than 99.8%. 
Then, the gate operation is performed, followed by an optional carrier pulse for analysis. Finally, 
we measure the probability p k of finding k ions in the \S) state by detecting light scattered on the 
5*1/2 ^ P\/2 dipole transition with a photomultiplier for 3 ms. The error in state detection due to 
spontaneous emission is estimated to be less than 0.15%. Each experimental cycle is synchronised 
with the frequency of the AC-power line and repeated 50-200 times. The laser beam performing 
the entangling operation is controlled by a double-pass acousto-optic modulator (AOM) which al- 
lows setting the frequency uo L and phase of the beam. By means of a variable gain amplifier, we 
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control the radio-frequency (r.f.) input power and hence the intensity profile of each laser pulse. 
To generate a bichromatic light field, the beam is passed through another AOM in single-pass con- 
figuration that is driven simultaneously by two r.f. signals with difference frequency 5/n. Phase 
coherence of the laser frequencies is maintained by phase-locking all r.f. sources to an ultra-stable 
quartz oscillator. We use 1.8 mW of light focused down to a spot size of 14 pm Gaussian beam 
waist illuminating both ions from an angle of 45° with equal intensity to achieve the Rabi frequen- 
cies Q/(2tt) ~ 110 kHz required for performing a gate operation with [u — 8)/ (2ir) = 20 kHz and 
r\ = 0.044. To make the bichromatic laser pulses independent of the phase £, the pulse is switched 
on and off by using pulse slopes of duration t t = 2 /as. 

Multiple application of the bichromatic pulse of duration r gate ideally maps the state \SS) to 
\SS) ^ \SS) + i\DD) ^ \DD) ^ \DD) + i\SS) ^ \SS) ^ ... (2) 

S v ' 

#1 

up to global phases. Maximally entangled states occur at instances r m = m ■ r gate (m — 1, 3, . . .). 
A similar mapping of product states onto Bell states and vice versa also occurs when starting 
from state \SD). In order to assess the fidelity of the gate operation, we adapt the strategy 
first applied in refs. 11 ' 13 consisting in measuring the fidelity of Bell states created by a single 
application of the gate to the state \SS) (see Fig. [2] (a)). The fidelity F = (^i|p eip |^i) = 
(pTsss + Pdddd)/^ + I m PDD ss> wmi me density matrix p exp describing the experimentally 
produced qubits' state, is inferred from measurements on a set of 42, 400 Bell states continuously 
produced within a measurement time of 35 minutes. Fluorescence measurements on 13, 000 Bell 
states reveal that p e s X s t ss + Pdd,dd = P2+P0 = 0.9965(4). The off-diagonal element p^D,ss 1S 
termined by measuring P(4>) = (cr^cr^ ) for different values of 0, where = a x cos 4> + cr y sin 0, 
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by applying (^-pulses to the remaining 29, 400 states and measuring p + p 2 — p\ to obtain the 
parity (cr^a^). The resulting parity oscillation P(4>) shown in Fig.[2](b) is fitted with a function 
Pfit{4>) — ^4sin(20 + O ) mat yields A = ^Ip^d ss\ = 0.990(1). Combining the two measure- 
ments, we obtain the fidelity F = 99.3(1)% for the Bell state 

A wealth of further information is obtained by studying the state dynamics under the action 
of the gate Hamiltonian. Starting from state \ SS), Fig.|3]depicts the time evolution of the state pop- 
ulations for pulse lengths equivalent up to 17 gate times. The ions are entangled and disentangled 
consecutively up to nine times, the populations closely following the predicted unitary evolution 
of the propagator (OQ) for ( = shown in Fig.[3]as a solid line. 

To study sources of gate imperfections we measured the fidelity of Bell states obtained after 
a pulse length r m for up to m = 21 gate operations. The sum of the populations p (t) + P2(t) does 
not return perfectly to one at times r m as shown in Fig. |4]but decreases by about 0.0022(1) per 
gate. This linear decrease could be explained by resonant spin flip processes caused by spectral 
components of the qubit laser that are far outside the laser's linewidth of 20 Hz (see Methods). 
The figure also shows the amplitude of parity oscillation scans at odd integer multiples of T gate 
similar to the one in Fig. [2] (b). The Gaussian shape of the amplitude decay is consistent with 
low-frequency noise of the magnetic field and the laser frequency as the source of imperfections 
(see Methods). 

The observed Bell state infidelity of 7 ■ 10 -3 suggests that the gate operation has an infidelity 
below the error threshold required by some models of fault-tolerant quantum computation 6-8 . How- 
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ever, further experimental advances will be needed before fault-tolerant computation will become 
a reality as the overhead implied by these models is considerable. Nevertheless, in addition to mak- 
ing the implementation of quantum algorithms with tens of entangling operations look realistic, the 
gate presented here also opens interesting perspectives for generating multi-particle entanglement 22 
by a single laser interacting with more than two qubits at once. For the generation of N-qubit GHZ 
states, there exist no constraints on the positioning of ions in the bichromatic beam that otherwise 
makes generation of GHZ states beyond N=6 so difficult in the case of hyperfine qubits 23 . While 
the bichromatic force lacks a strong spatial modulation that would enable tailoring of the gate 
interaction by choosing particular ion spacings 24 25 , more complex multi-qubit interactions could 
be engineered by interleaving entangling laser pulses addressing all qubits with a focussed laser 
inducing phase shifts in single qubits. Akin to nuclear magnetic resonance techniques, this method 
should allow for refocussing of unwanted qubit-qubit interactions 26 and open the door to a wide 
variety of entangling multi-qubit interactions. 
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Methods 

AC-Stark-shift compensation The red- and the blue-detuned frequency components u± of the 
Dichromatic light field cause dynamic (ac-) Stark shifts by non-resonant excitation on the carrier 
and the first-order sidebands that exactly cancel each other if the corresponding laser intensities 
I± are equal. The remaining ac-Stark shift due to other Zeeman transitions and far-detuned dipole 
transitions amounts to 7 kHz for a gate time T gate = 50 /is. These shifts could be compensated by 
using an additional far-detuned light field 27 or by properly setting the intensity ratio We 
utilize the latter technique which makes the coupling strengths Qss^+dd oc 2 v / / + /_, £1 S d^>ds oc 
I + + J_ slightly unequal. However, the error is insignificant as ^Isd^ds/^ss^dd — 1=4- 1CT 3 
in our experiments. 

Sources of gate infidelity Spin flips induced by incoherent off -resonant light of the bichromatic 
laser field reduce the gate fidelity. A beat frequency measurement between the gate laser and a 
similar independent laser system that was spectrally filtered indicates that a fraction 7 of about 
2 ■ 1CT 7 of the total laser power is contained in a 20 kHz bandwidth B around the carrier transition 
when the laser is tuned close to a motional sideband. A simple model predicts spin flips to cause 
a gate error with probability pfu p = [ti^\v — S\)/ (2rfB). This would correspond to a probability 
PfUp = 8 • 10~ 4 whereas the measured state populations shown in Fig. |4] would be consistent with 
Pfii P = 2 • 10~ 3 . Spin flip errors could be further reduced by two orders of magnitude by spectrally 
filtering the laser light and increasing the trap frequency u to above 2 MHz where noise caused by 
the laser frequency stabilization is much reduced. Variations in the coupling strength induced 
by low-frequency laser intensity noise and thermally occupied radial modes can be independently 
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measured by recording the amplitude decay of carrier oscillations. For 5Q / Q, we find a variation of 
7 • 10~ 3 . The Gaussian decay of the parity contrast shown in Fig. © is attributed to low frequency 
noise randomly shifting the laser frequency with respect to the atomic transition frequency lu and 
having a Gaussian distribution with full-width at half maximum of 5co/(2tt) = 180 Hz. This value 
is consistent with Ramsey measurements on a single ion predicting 5u/(2tc) = 160 Hz. For a 
single gate, the frequency uncertainty gives rise to a fidelity loss of 0.1%. 

A bichromatic force with time-dependent f2(t) acting on ions prepared in an eigenstate of S y 
creates coherent states a(t) following trajectories in phase space that generally do not close 19 ' 28 . 
For the short rise times used in our experiments, this effect can be made negligibly small by slightly 
increasing the gate time. 
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Figure 1: Gate mechanism, (a) A bichromatic laser field with frequencies lu + , lu_ satisfying 2u 
= uo + + uo_ is tuned close to the upper and lower motional sideband of the qubit transition. The 
field couples the qubit states \SS) <-> \DD) via the four interfering paths shown in the figure, 
n denoting the vibrational quantum number of the axial COM mode. Similar processes couple 
the states \SD) <-> \DS). (b) The qubits are encoded in the ground state Si/ 2 (m = 1/2) and 
the metastable state D 5 / 2 (m = 3/2) of 40 Ca + ions and manipulated by a narrow bandwidth laser 
emitting at a wavelength of 729 nm. 




|SS,n> 
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Figure 2: High-fidelity gate operation, (a) State evolution induced by a M0lmer-S0rensen bichro- 
matic pulse of duration r. The Rabi frequency £l(t) is smoothly switched on and off within 2 jis 
and adjusted such that a maximally entangled state is created at r gate = 50 fxs. The dashed lines 
are calculated for h CO m = 0.05 from the propagator (OQ), neglecting pulse shaping and non-resonant 
carrier excitation. The solid lines are obtained from numerically solving the Schrodinger equa- 
tion for time-dependent Cl(t) and imbalanced Rabi frequencies = 1.094 (see Methods), 
(b) A (|)^ analysis pulse applied to both ions prepared in ^ 1 gives rise to a parity oscillation 
P{4>) = sin(20) as a function of <p. A fit with a function Pf it = A sin(20 + <f> ) yields the parity 
oscillation amplitude A = 0.990(1) and (fio/ir = —1.253(1). The precise value of the phase 4> is 
without significance. It arises from phase-locking the frequencies lu , lu + , io- and could have been 
experimentally adjusted to zero. 
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Figure 3: Entanglement and disentanglement dynamics of the M0lmer-S0rensen interaction. 

Starting from state \SS) for a detuning of the bichromatic laser from the sidebands set to S — v = 
—20 kHz, the figure shows the time evolution of the populations po, pi, and p 2 denoted by the 
symbols (•), (o), and (A), respectively. The length of the pulse is equivalent to the application 
of up to 17 gate operations. Maximally entangled states are created whenever po(r) and P2{t) 
coincide and p 1 (r) vanishes. 
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Figure 4: Multiple gate operations. Gate imperfections as a function of the bichromatic pulse 
length r m = m ■ T gate given in equivalent number of gate operations m. The upper curve shows 
a linear decrease of the state populations p + p 2 with a slope of 0.0022(1). The lower curve 
the magnitude of the coherence 2pDD j ss measured by detecting parity oscillations and fitted by 
a Gaussian decay function that accounts for low-frequency noise of the laser frequency and the 
magnetic field. Combining both measurements yields the Bell state fidelity F m shown as the middle 
trace. For m = 21, the fidelity is still F 2 i = 80(1)%. Similar results are achieved when replacing 
the entangling pulse of length r m by m amplitude-shaped pulses each of which is realizing an 
entangling gate operation. 
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